Strong Scalar QED in Inhomogeneous Electromagnetic Fields 
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Strong QED has attracted attention recently partly because many astrophysical phenomena have 
been observed to involve electromagnetic fields beyond the critical strength for electron-positron pair 
production and partly because terrestrial experiments will generate electromagnetic fields above or 
near the critical strength in the near future. In this talk we critically review QED phenomena 
involving strong external electromagnetic fields. Strong QED is characterized by vacuum polariza- 
tion due to quantum fluctuations and pair production due to the vacuum instability. A canonical 
method is elaborated for pair production at zero or finite temperature by inhomogeneous electric 
fields. An algorithm is advanced to calculate pair production rate for electric fields acting for finite 
periods of time or localized in space or oscillating electric fields. Finally, strong QED is discussed 
in astrophysics, in particular, strange stars. 
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I. INTRODUCTION 

Recently strong QED (quantum electrodynamics) has 
attracted much attention not only from theoretical in- 
terest but also from astrophysical observations and ter- 
restrial experimental tests in the near future. From a 
theoretical view point, calculating the full nonperturba- 
tive effective action under the influence of strong exter- 
nal electromagnetic fields, in particular, inhomogeneous 
fields, is still a challenging task (for a recent review and 
references, see Ref. jlj and also Ref. [5|). From an ex- 
perimental view point, in the near future electromagnetic 
fields from X-ray free electron lasers from LCLS (Linac 
Coherent Light Source) at SLAC ^ and TESLA (TeV 
Energy Superconducting Linear Accelerator) at DESY 
Q may attain a strength almost comparable to the criti- 
cal value for electron-positron pair production, which will 
directly test strong QED 5] . Interestingly, astrophysical 
sources have been predicted and observed that can have 
electromagnetic fields greater than the critical strength. 
Neutron stars have magnetic fields ranging from 10^ G to 
10^^ G and more than one-tenth of them have magnetic 
fields stronger than 10^'' G, the so-called magnetars (for 
a review and references, see Ref. @), at least one order 
greater than the critical strength. Another interesting 
astrophysical objects with a ultra-strong electromagnetic 
field are strange quark stars, hypothetical objects, which 
may have electric fields with one or two order higher than 
the critical strength 0,i] (see also Ref. [9,] for review and 
references) . 

Vacuum fluctuations due to a strong external electro- 
magnetic field contribute nonlinear terms to the classi- 
cal Maxwell theory and the electromagnetic theory thus 



becomes highly nonlinear. Physics in strong electromag- 
netic fields drastically differs from the Maxwell theory 
[loL |Tl[ . The cyclotron energy of an electron in a strong 
magnetic field can be greater than the rest mass energy 
of electron, the equivalent value leading to the critical 
strength of magnetic field Be = m'^c^/eh (4.4 x 10^'^ G). 
Similarly, in a strong electric field, virtual pairs of elec- 
trons and positrons can gain energy comparable to or 
greater than the rest mass energy of electron or positron. 
The electric field whose potential energy across the 
Compton wavelength is the rest mass energy of electron is 



the critical value En 



,2^3 



/eh (2.2 X IQis V/cm). For 
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magnetic fields greater than the critical value, nonlin- 
ear contributions to the Maxwell term make the vacuum 
polarized by quantum fluctuations and the vacuum po- 
larization causes nonlinear effects such as birefringence 
(propagation of photons in the magnetic vacuum) , which 
plays an important role in the physics of magnetars [IQ] . 
For strong electric fields the vacuum decays due to an 
imaginary part of the effective action and thus leads 
to Schwinger pair production [l^. Strange stars can 
emit electron-positron pairs more efficiently than pho- 
tons [H, [3, On the other hand, in the standard 
QED with the minimal interaction magnetic flelds are 
stable up to B = 10'^^ G due to the instability from 
the self-interaction of an electron and up to the range 
B = 10^^ — 10^^ G due to the instability from magnetic 
monopole production at the string or Planck scale p^ . 
However, the Pauli interaction may open a window for 
pair production by a far weaker inhomogeneous r nag netic 
field and would have astrophysical applications [17| . 

QED describes the interaction between charged parti- 
cles and photons. The success of QED is based on the 
perturbation theory in the weak-field limit. However, 
QED has not been completely understood yet in the op- 
posite case of strong electromagnetic fields partly because 
the full nonperturbative QED action is not known ex- 
cept for some exactly solved cases [l|, Q • Historically, the 



2 



effective action of an electron in a constant electromag- 
netic field was obtained by Heisenberg and Euler [l^ , and 
also by Weisskopf [l^. Using the proper time method, 
Schwinger found the one-loop effective action for a spin- 
1/2 fermion with charge q and mass m in a constant 
electromagnetic field |12l | 



-off 



1 



87r2 



ds- 



2 Ke cosh{qsX) 



1 



'-{gsfT . (1) 



Imcosh(gsX) 

where 

X = [2{J^ + ig)]^^^ = Xr + iX, 
Here, J- is the negative of the Maxwell term, — /^MaxwcU, 
1„ ^„., 1, 



(2) 



and g is another Lorcntz invariant tensor 



E B, 



(3) 



(4) 



where F^"' = et^^^l^Fafs is the dual field tensor. The one- 
loop effective action was also obtained in Ref. [20| . 

The one-loop effective action has two important as- 
pects. First, in the weak-field limit the nonlinear contri- 
bution to the real part 



Re/:(i 



(5) 



makes the vacuum polarized by quantum fluctuations. In 
a pure strong magnetic field the leading term becomes 



£(1) 



{qBf (2qB^ 



(6) 



The ratio r^^V-CMaxwdi = -{q^ /l2T:^)\n{qB /m^) is the 
one-loop QED /3-function related with the renormaliza- 
tion group [1]. Second, in spinor QED in a pure strong 
electric field pairs are produced at the rate per unit time 
and unit volume 



..fcrr 



' n=l 



(7) 



and in scalar QED pairs at the rate 



vacuum solution of the field equation and then calculate 
the number of pairs produced by the field. The tunnel- 
ing interpretation is that charged pairs in the Dirac sea 
can tunnel quantum mechanically through the potential 
barrier lowered by a uniform electric field. Pair produc- 
tion is an efficient mechanism for energy extraction from 
objects with strong electromagnetic fields. 

In this talk, within the framework of canonical quan- 
tum field theory we critically review the Schwinger mech- 
anism at zero or finite temperature in inhomogeneous 
electric fields motivated by terrestrial experiments or as- 
trophysics. Exact solutions of the Klein-Gordon or Dirac 
equation minimally coupled to inhomogeneous electro- 
magnetic fields are known only for a few models. The 
Sauter-type electric field that extends for a finite region 
or lasts for a finite period is the most well-known model 
[2l[. However, one cannot find, in general, solutions for 
arbitrary electromagnetic fields, so he or she has to em- 
ploy some approximation schemes. It is known that pair 
production by a constant electric field that extend over a 
finite region has a finite size effect and differs from that 
by a constant field (2^ . Also the pair production rate 
by a Sauter-type electric field obtained by the worldline 
instanton method depends on the characteristic scale in 
a nontrivial way |23l . [2^ . Applying the phase-integral 
method [1^ to find the WKB instanton action for the 
field equation with an electric and/or magnetic field in a 
fixed direction, the pair production rate is obtained for 
the Sauter-type electric field cither in space or time with 
or without a constant magnetic field [26,, 2lL, 2M- Fur- 
ther, a perturbative method is advanced to calculate the 
WKB instanton action for pair production by any ana- 
lytical electric field and is then applied to strange stars to 
calculate the production rate of electron-positron pairs. 
The thermal effect on the pair production is also studied 

The organization of this talk is as follows. In Sec. II, 
we critically review the Schwinger mechanism and then 
thermal effects on pair production. In Sees. HI and IV, 
we apply the WKB instanton action method to inhomo- 
geneous electric fields that act on for a finite period or 
extend for a finite region or oscillate. In Sec. IV, we 
apply strong QED to calculate the pair production rate 
from strange stars. 



II. CANONICAL METHOD FOR PAIR 
PRODUCTION 



yboson ^ 2Im/:|?^'^°" 



^ ' n=l 



(-1) 



.(8) 



Here, the factor of 2 is the spin multiplicity. The scat- 
tering amplitude of the ingoing vacuum to the outgo- 
ing vacuum decays according to Eqs. ([7]) or ([8]), leading 
to Schwinger pair production. One way to understand 
Schwinger pair production is to compute the imaginary 
part of the effective action. Another way is to find the 



In real physical systems electric fields are either con- 
fined to finite regions or turned on for finite periods of 
time. For such electric fields it is a nontrivial task to cal- 
culate the pair production rate. Instead of applying the 
proper time method or path integral method, we employ 
an approximation scheme such as the WKB approxima- 
tion and phase integral in canonical quantum field theory. 
For the sake of convenience we consider only scalar QED, 
but the formalism here can be directly applied to spinor 
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QED [27|,l28 



The inverse Bogoliubov transformations are 



A. Schwinger Pair Production 

The Klein-Gordon equation for a charged boson with 
q (q > 0) and m takes the form (in units with h = c = 1 
and with metric signature (+, — , — , — )) 

[ri'"'id^ + iqAf,)id,+iqA,)+m^]^K,t)^0. (9) 

Hereafter we further restrict our study to time-dependent 
electric fields along the z direction with gauge potentials 
of the form Az{t) = —Eog{t) for any analytic function 
g{t). Then the Fourier mode, $(x, — e''^'^(f>]f^, satisfies 

.2 + ki + {k, + qEog{t)f] ^k(0 = 0. (10) 



9*2 



The solution can be used to quantize the position opera- 
tors as 



(11) 



(12) 



kit) = ^i,{t)ai,{t) + ^l{t)bl{t), 
and the momentum operators as 

K{t) = ^*(t)fet(t) + ^k(t)ak(i). 



On the other hand, in quantum mechanics, Eq. (jlO[) is 
a one-dimensional scattering problem with inverted po- 
tential. The positive (asymptotic) solution (^k,in at one 
asymptotic region t — —oo defines the (asymptotic) in- 
going vacuum and another positive (asymptotic) solution 
V'k.out at the other region t = oo defines the (asymptotic) 
outgoing vacuum. As an incident solution from t — oo 
is partially transmitted over the barrier to t = — oo and 
partially reflected by the barrier back to * = oo, the in- 
going solution is related with the outgoing solution as 



V'k,!!! = A^kV'k.out + fk'/'k out- 



(13) 



That is, the ingoing positive frequency solution is mixed 
both with the outgoing positive solution and with the 
outgoing negative solution, which is the origin of particle 
production by an external field [13, HH . As the Wron- 
skian 

^mv^.{t)-^^{t)^i{t) = i, (14) 

is constant, the coefficients satisfy the relation 

W?-W? = l. (15) 

In fact, the annihilation and creation operators at two 
asymptotic regions are related through Bogoliubov trans- 
formations 



ak,in - Mk^k.out - '^k^k.out' 



Ok, out = Mkak,in + '^k&kjr 
&k,out = Aik^k.in + i^kfllcii: 



(17) 



Therefore, the outgoin g va cuum contains the ingoing par- 
ticles/antiparticles as pol IsH 

(O,out| ^aj^_i^ak,i„|0,out) = ^ |i^kP, 

k k 

(0,out|^S[i„Sk,i„|0,out) = ^li^kP, (18) 

k k 

and, conversely, the ingoing vacuum evolves into outgoing 
particle/antiparticle states as 

(O,in|^aj^_^^jak,out|0,in) = ^ ji'kP, 

k k 

(0'in|^S],„Jk,out|0,in) = J^li^kp. (19) 

k k 

B. Hamiltonian Approach 

Pair production by time-dependent electric fields can 
also be described by the Hamiltonian formalism. The 
Hamiltonian formalism is particulary appropriate for 
studying thermal effects because the density operator 
should satisfy the Liouville-von Neumann equation with 
respect to the Hamiltonian itself. The complex scalar 
field has the Hamiltonian 

(20) 



H = j [rfk] TT^TTk + CJk(^)<^k'/'k 

where [c?k] = (P'k/{2T:Y and 

ujl{t) = + ki + [k, + qEog{t)f. 
The field operators are quantized as 

0*(i,x) = J[dk][Mt)k{t) + ipUt)aiit) 
and the momentum operators as 

n{t, ^) = j [dk] [^l{t)al{t) + (^k(t)6k(t) 
^*(i,x) = /[dk]U;;(i)S[(i) + (^k(i)ak(t) 



(21) 



.-ik-: 



„ik. -x 



(23) 



In the Hamiltonian approach both quantum states and 
the density operator can be found simultaneously us- 
ing the operators that satisfy the Liouville-von Neumann 
equation [H, fssj 



Mk^k, 



_ * -t 

out ^k^k,ouf 



(16) 



(24) 
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In fact, there are the time-dependent annihilation and 
creation operators satisfying Eq. ([M]) for particles 



(25) 



and for antiparticles 



(26) 



Then the time-dependent vacuum, an exact state of the 
time-dependent Schrodinger equation, is given by 

a^{t)\0;t) ^bk{t)\0;t) ^0 (for any k), (27) 

and multi-particle and antiparticle states by 



I'^ki 



;"-k2 



47 w ^irw 



/ni! 



\0;t). (28) 



C. Pair Production at Finite Temperature 

The finite temperature QED effective action was cal- 
culated in a constant magnetic field ^33| , a constant elec- 
tromagnetic field |35] , and at finite density [36]. The 
Schwinger proper-time method was used to derive the ef- 
fective action in a constant electromagnetic field, which 
exhibits the Schwinger mechanism at high temperature 
[37| . However, depending on the formalism employed to 
calculate the effective action, pairs are either produced 
[aSllil or not produced The QED effective action 
from the imaginary-time formalism has nonzero imagi- 
nary part at two- loop (41i] . In this paper we follow the 
real-time formalism in Ref. [29| to obtain pair production 
at finite temperature at one-loop. 

As ak(t),ak ^^'^ ^k(i),^k satisfy Eg. (l24t . the density 
operator for particles can be found as [33| 



PakW = ^exp 

^k 



-/3<(4W«k(t) + i)], (29) 



with (3 = l/(fcT) is the inverse temperature, and for an- 
tiparticles as 



1 



PbAt) = ^ exp -/3<( 6^ (t)fek(O + 7^ 



1 



(30) 



Then the pair production rate for each spin component 
from an initial ensemble by an electric field acting for a 
finite period of time is given by 29] 



n^iE,T) = Tr(pif7V°f)-/^ 

= \ME)\'{VtiT) + l), 
where f^^ is the Bose-Einstein distribution 

1 



1 



(31) 



(32) 



Here, the Bogoliubov coefficients are 



t'k 



(oo) = z(v.^(oo)^r-vi^(oo)v.r). (33) 



When there is a uniform magnetic field B in addition to 
E(t), each mode of the scalar field obeys the equation 



ipnkit)+u;lk{t)ip^kit) = 0, (34) 



where 



4fc = {k, + qEog{t)f + qB{2n + 1) + (35) 



III. TIME-DEPENDENT ELECTRIC FIELDS 

Pair production by localized electric fields in time 
or space significantly differs from that by the con- 
stant electric field due to a duration or a size effect 
m m [13, m, 113. in this section we exploit an analyt- 
ical method to calculate the Schwinger pair production 
rate by an electric field acting for a finite period of time 
or an oscillating electric field. This case is characterized 
by a homogeneous time-dependent electric field E{t) with 
the maximum strength Eq and the time scale T defined 
as 



T : 



1 



E{t)dt. 



(36) 



The pair production rate is determined by two dimen- 
sionless parameters 



qEo 



qEoT' 



(37) 



Pair production is allowed for any e but is strongly sup- 
pressed for e 1. 

The main result of Ref. is that in the weak-field 
limit {E < Ec) the mean number of boson pairs for each 
mode k per unit time and unit volume 

AAk = e--^-- (38) 

is determined by the WKB instanton action of Eq. (fTU]) 

Sw^ij) [k, + qEog{t) f + m2 + k^di, (39) 

where the integral is taken outside the contour in the 
complex plane of time. In Ref. [T^ the production rate 
in scalar QED is defined as twice the imaginary part of 
the effective action 



Wk = 21n(l + 6"'^'') =2^ 



(-1) 



n+l 



-nSk 



(40) 



The production rate (j40p is even valid for very strong 
electric fields [E > E^) or (5k < 1). 
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A few comments are in order. Taking into account spin 
statistics, the mean number of boson pairs with spin s is 
given by [il ] 



^^boson 



1 = (1 + e 



-5£ 



\2s+l 



1, (41) 



while the mean number of fermion pairs with spin s 



■fermion 



1 - e-""^ 



1 - (1 - 



n Doson/fcrmion ■ i ■ i i ■ ■ j_ r j_i r 

where w^^ is twice the imaginary part oi the ef- 

fective action for bosons and fermions, respectively, and 
2s+l is the spin multiplicity. The mean numbers without 
the spin multiplicity is also obtained in Ref. [i^. The 
mean numbers (j4ip and (|42p hold even for very strong 
electric fields. The first term in Eq. ((4T|) is the ampli- 
fication factor for boson production. In the weak-field 
limit (iSk > 1), the mean number of boson or fermion 
pairs is approximately given by 



boson/fcrmion 



(2s + l)e' 



n/fen 



(43) 



Another interesting point is the relation between the 
WKB instanton actions for bosons and fermions. For 
a Sauter electric field it is shown in Ref. 28] that the 
WKB instanton action for fermions together with the 
next-to-leading order is nothing but the WKB instanton 
action for bosons. Such a relation plausibly holds for any 
inhomogeneous electric field since the one-loop effective 
action for fermions 



S 



'fermion 



= -ilndet(iD - m), D = 7^9^ + iqAf,, (44) 



has the imaginary part, which is approximately for 
bosons 



Im(5 



fermion^ ^ - - [In det(iD - m) + In det(iD - m)*] 



--lndet(i:)2 +to2). 



(45) 



Therefore, the WKB instanton action works for both 
scalar and spinor QED. 

Now we develop an algorithm to compute the instanton 
action systematically. Introducing another variable 



C = 5(0, 



(46) 



we rewrite Eq. ([55]) as 



Sk = iiqEo) S Ml 



and expand the square root in an inverse power series 



n=0 



and the function l/g\C,) in a power series 

^ oo 

7(0 



(49) 



Then the sum of negative residues of simple poles leads 
to the WKB instanton action 



5k = 27r(g£'o) ^ C„+2l'„ 



(50) 



n=0 



The first few terms of C„ are 

Co = 1, Ci = ai, C2 = 



Q!2 



C3 = ^^^^-y ^^1) 



where 



ai 



-, 012 



qE^ ' ' - {qE^Y 



(52) 



The coefficients I?„ are determined by the profile of g{t). 
For specific models, we consider a Sauter-type electric 
field E{t) — Eq sech^(i/T) and an oscillating electric 
field E{t) = Eo cos{t/T). 



A. Sauter- Type Electric Field 



The gauge potential for E{t) — Eq sech^ (t/T) in the 
direction is given by the Sauter potential 



Mt) 



^So^tanhlJ; 



(53) 



With the change of variable C, = g{t) — rtanh(i/T), we 
have the power series 



^ ' ^ W n=0 



(54) 



and find 2?2n = l/T^". This means that the WKB in- 
stanton action is 



00 ^ 

S. = 2.iqEoT^)j:^. 

n=Q 

In terms of the scaled variables and parameters 



A 



k 



^ Z ^ 2nqEoT^ 



qEoT m 
the leading terms of the WKB instanton action are 

-1 e\l + K^) 

-2^2 8 



(55) 



, (56) 



5k = Ze^il + iY 



(57) 



In fact, the sum ((55|) can be done exactly as [25 
Z 



V/(1+A)2+ £2(1 + ^2) 



+ \/(l-A)2+e2(i + K2)_2 



(58) 



The instanton method ^281 gives much closer result to the 
exact one [21] than the worldline instantons [23 |. 



6 



B. Oscillating Electric Field 



The oscillating electric field E{t) — Eq cos{t/T) has 
many physical applications such as laser fields [ij, [H, 
l46l | . The electric field from oscillating plasma due to pair 
production is approximately given by E{t) — E cos(t/T) , 
where E varies slowly during the oscillation period |47l . 
m, m, The gauge potential C = git) = Tsm{t/T) 
leads to the expansion 



1 



9'it) 



= 1 



1- ^ 

^ T2 



2r2 ^ 8T4 ^ lere 



.(59) 



Repeating the procedure in Sec. Ill A, we obtain the 
leading terms of the WKB instanton action 



5k = Ze^{l 



ri 

2 



A2 _ e^jl + n') 
4 " 16 



(60) 



IV. INHOMOGENEOUS ELECTRIC FIELD 

For an inhomogeneous electric field localized in the 
z direction, we may choose a Coulomb gauge Ao{z) = 
—EQh{z), which leads to E{z) = Egh'^z). Then the 
Klein-Gordon equation has the Fourier mode solution, 



<I>(x, t) 



„ikj_ ■x_L —iut. 



>k^ , given by 



The difference from the case of time-dependent electric 
fields is that the overall sign changes, and ai and a2 are 
now replaced by ai = oj/iqEo) and 012 — —a2- 

For specific models, we consider first a Sauter-type 
electric field E{z) = Eq secYi' [z/L) and then the electric 
field from strange stars in the next section. The Coulomb 
gauge is given by the Sauter potential 



Aq{z) = -EQLiSinh(^^. 



(67) 



With Q — Ltanh(2/i) and the replacement of ai by ai 
and a2 by a2 = —02, the leading terms of the WKB 
instanton action are 



5k = Zf{l + K^) 
where 



1 A2 £2(1 ^ ^2) 

- H \ 

2 2 8 



(68) 



A = 



, KEE— , Z = 2nqEoL^ ^ —r. (69) 



qEoL m 
The exact sum of (1551) is known [2 



5k, 



(1 + A)2+ £2(1 + ^2) 



-V(l-A)2+ £-2(1 + ^2) 



(70) 



Q2 

-— + m2 + ki - (c^ + qEog{z))'' (^k^ (z) = 0. (61) 



The characteristic length scale L is defined as 



1 /"^ 

^-C'O J-00 



(62) 



As for time-dependent electric fields, two dimensionless 
parameters 



qEo 



qEoL' 



(63) 



determine the pair production rate. Remarkably the pair 
production rate is again given by the WKB instanton 
action (H] 

S^ = -i<f ^{uj + qEoh{z) f - (to2 + k2_)dz, (64) 



where the integral is taken outside the contour in the 
complex plane of space. We again introduce the variable 



C = /i(z), 
and rewrite Eq. ([64|) as 



5k = -i 



2 m2+k2^ C^C 
" {qEoCr 5'(C) 



(65) 



(66) 



V. PAIR PRODUCTION FROM STRANGE 
QUARK STARS 

A source of the most strong electric fields beyond the 
critical strength is strange quark stars. A quark star, a 
hypothetical astrophysical object, could be formed from 
a hadron-quarkphase transition at high densities and/or 
temperatures jsil . [s^ . [53| . Chemical equilibrium of u, d 
and s quark and the charge neutrality of strange stars 
requires a net amount of electrons that are free to move 
the surface but bounded by the electric attraction from 
the positive core, thus forming an electrosphere of several 
hundred fermis. Then, the electrosphere of strange stars 
can generate an extremely strong electric field as strong 
as 5 X 10^^ V/cm, two order greater than the critical 
strength, and leads to efficient production of electron- 
positron pairs [13., J^, [l^ . 

The static potential from the Thomas- Fermi model of 
the electron distribution at temperature T [7| 



(y3-y/) + 7r2T2(y-l/,) , (z<0). 



d^V 


4a ■ 


dz^ 


3^ . 


d^V 


4a ■ 


dz^ 


3^ . 



> 0), 



(71) 



where z is the coordinate normal to the quark surface, 
a is the fine structure constant, and Vq/Sir'^ is the quark 
charge density inside the quark matter. The boundary 
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condition is V(-oo) = Vq, V{oo) = 0, and y(0) = 3K,/4. 
The Coulomb gauge potential is found [l^ 



sinh 

and the electric field is 



(72) 



E{z,T) 



87r3^2 cosh[2y^T(z + ^o)] 



3 sinh2[2y^T(z + zo)]' 
whose characteristic scales are 



So 



L 



(73) 



(74) 



The higher (lower) the temperature is, the greater 
(smaller) is the maximum strength and the narrower 
(wider) is the width of the electric field. 
With 



c 



L 



sinh 



(75) 



the WKB instanton action becomes 




Thus, in terms of the scaled variables and parameters in 
Sec. IV, the leading terms of the WKB instanton action 
are 



5k = 27r{qEoL'') 

= zx\i 



1 u) + 



qEoL 4 qEoL {qEoLf 
e2(l + «2) 



(77) 



Then, the mean number (|42p of electron-positron pairs 
per unit time and volume 



^^fcrmion 



(78) 



is the spectrum of emitted pairs. As Z\ — ^Jiir /a x 
{uj/T), hot strange stars produce more pairs of electrons 
and positrons than cold ones, confirming the numerical 
result of Ref. 



electromagnetic fields are known only for a few cases, for 
general electromagnetic fields, however, one has to rely 
on some approximation schemes. Inhomogeneous electric 
fields result in a finite size or duration effect and differs 
from that by a constant field [l^, [2^ [13] ■ We applied 
the phase-integral method to find the WKB instanton 
action for the mode equations in inhomogeneous electro- 
magnetic fields and then calculated the pair production 
rate by a Sauter-type electric field either in space or time 
[H, [23, [ill and an oscillating electric field. We also stud- 
ied the thermal effect on pair production by an electric 
field that acts for a finite period of time [29]. Finally, 
we applied the WKB instanton action method to strange 
stars to calculate the electron-positron pair production 
rate. 

The issues not treated in this talk are the effective 
action and the back reaction of QED at zero or finite 
temperature. It is a complicated task to obtain the effec- 
tive action in inhomogeneous electromagnetic fields. In 
canonical quantum field theory, we may follow Ref. [s^ ]. 
according to which the effective action is related with the 
scattering amplitude as 



(O,out]0,i 



(79) 



Thus the effective action requires a complete knowledge 
of evolution of the ingoing vacuum to the outgoing vac- 
uum, which may follow from the vacuum wave functional 
from Eqs. (P5)) and (|26p for each mode. Another impor- 
tant issue is the QED back reaction problem, which is 
described by, for instance, the scalar QED action of the 
form 

C = r{d^ + ^qA^fcf>-^F^,F''r (80) 

The back reaction cannot be neglected for strong electric 
fields because the additional electric field produced by 
pairs is comparable to the applied field. In fact, positive 
(negative) charges of produced pairs move in the same 
(opposite) direction of the applied electric field, so the 
current due to pairs induces an electric field in the oppo- 
site direction of the applied field and overshoots it until 
the process is reversed, which leads to the famous plasma 
oscillation [13, [H, US', 50^ • These issues will addressed in 
a future publication ^55.] . 
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